Spin-density induced by electromagnetic wave in two-dimensional electron gas with 
both Rashba and Dresselhaus spin-orbit couplings 
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We consider the magnetic response of a two-dimensional electron gas (2DEG) with both Rashba 
and Dresselhaus spin-orbit coupling to a microwave excitation. We generalize the results of Ref. [lj], 
where pure Rashba coupling was studied. We observe that the microwave with the in-plane electric 
field and the out-of-plane magnetic field creates an out-of-plane spin polarization. The effect is more 
prominent in clean systems with resolved spin-orbit-split subbands. Considered as response to the 
microwave magnetic field, the spin-orbit contribution to the magnetization far exceeds the usual 
Zeeman contribution in the clean limit. The effect vanishes when the Rashba and the Dresselhaus 
couplings have equal strength. 
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PACS numbers: 72.25.Rb,85.75.-d 



The spin-orbit effects in semiconductors have been 
studied for a long time 0, G3> El- The discussion was 
revived in relation to the spin-Hall effect Q in hope of 
applying spin-orbit related effects to spintronics Q. Ini- 
tially the effect was considered for conductors, e.g., for 
a 2DEG. Although the attention of the community has 
now mostly switched to the quantum spin-Hall effect in 
" spin-Hall insulators" [H, [t| , we consider here the spin- 
Hall related effects in a 2DEG with both Rashba and 
Dresselhaus spin-orbit couplings. 

The promise of the spin-Hall effect is in the possibil- 
ity of generation of non-equilibrium spin polarization by 
means of a DC electric field. However, after some discus- 
sions it was concluded that in two-dimensional electron 
gases with the Rashba and the Dresselhaus spin-orbit in- 
teractions the spin-Hall effect vanishes for constant and 
homogeneous electric field 

mind 

llal. At finite fre- 
quencies the spin-Hall effect is non-zero [111 ]. 

While for the homogeneous spin-Hall effect (uniform 
applied electric field) the out-of-plane spin polarization 
is expected to accumulate only at the edges of the sam- 
ple llj , there is an alternative possibility that we explore 
here: To create an out-of-plane inhomogeneous spin den- 
sity in the bulk in response to a spatially-modulated field. 
Such bulk accumulation is free from the uncertainties as- 
sociated with the charge and spin transport near the sam- 
ple boundaries, and thus may provide an unambiguous 
method to detect the spin-Hall effect. 

Alternatively, one can consider our results as providing 
the spin response to a long wave (spatially homogeneous), 
out-of-plane, oscillating magnetic field. According to the 
Faraday induction law such a field creates spatially in- 
homogeneous, in-plane electric field which in turn is re- 
sponsible for the effect. 

In this paper we generalize the results of Ref. [ij where 
the magnetic response of a 2DEG with pure Rashba spin- 



orbit coupling to microwaves was studied. Here we con- 
sider a more general situation where both couplings are 
present. 

The spin-orbit coupling reads 

Hso = <XR(-Px<r y + Py<Jx) + a D (p x cr x - Py<T y ) , (1) 

where aR and «d are the strengths of the Rashba and 
Dresselhaus spin-orbit couplings respectively. It is con- 
venient to perform a 7r/4 rotation in both the momen- 
tum and the spin spaces. That is p' x = Px+ J! v , p' 



and ai 



V2 



Px+Py 

" i T Then the 



spin-orbit coupling reads 



V2 



H S o = -{ctR + ar>)p' x a' y + (or - a^)p' y a' x . (2) 

In what follows we work in the rotated basis and omit 
the primes. 

Introducing the angle cj> p via p = |p|(cos</> p ,sin</> p ) 
we obtain the energies of the two sub-bands given by 
' — J^pr ± a^((/)p)\p\, where m* is the electron 
y/(a^ + al,)(l +T Q cos20 p ), and 



e^(p) 

band mass, a c f > (4> p ) 



2ttR«D 

di+a 



R 



g-. For purely Rashba (Dresselhaus) coupling 



T a — 0, while in the case of equal coupling strengths 
T a = 1. 

We consider a linearly polarized in-plane mi- 
crowave field A = Aoexp(iqr — iClt), where Aq = 

(cos a, sin a, 0) and q = q(sma, — cos a, 0). The signs 
are chosen so that for positive Aq and q the vectors 
q, Aq, e z form a right-handed basis. As usual E = 
(ifl/c)A and B = iq x A. 

Kinetic equation. We follow the route of Refs. P, [111, 
03] and use the standard linear response Keldysh tech- 
nique to determine the dynamics of the charge and spin 
densities. For introduction into the Keldysh technique 
see Ref. (H)]. In the dirty limit (to be defined below) 



2 



this technique leads to diffusion equations for the spin 
and charge densities 11, [l6j]- We concentrate here on 
the clean limit although our results are valid also in the 
dirty limit. 

We consider only the s-wave disorder scattering, that is 
Vdisordcr = J2k u ^i r ~ r k) where Tfc are random locations 
with the average density n- lmp . We employ the linear 
response, H = Hq + Hi, with 



H = 



P 



2m* 



VP 



^disorder ; 

where r) = [(-(a R + a D )cr y , (a R - a B )a x ], and 



(3) 



where v = 

m* 



2c 



{v,A}_ 



1 



-gp B B(7 , 



V- 



The zeroth order in A Green's functions, Go, reflect the 
standard disorder broadening. We introduce the inverse 



momentum relaxation time r 1 = 27 



''imp L 



n>imp is the density of impurities and v — m*/(2Trh ) is 
the electronic density of states per spin (strictly speak- 
ing v is the density of states in absence of the spin-orbit 
coupling). We obtain 



1 1 T)p 



(~iR-\ 



1 rjp 



2a<f,\p\ 



G R ~ , (4) 



. In equilib- 
tanh^|^. 



where G^fato) = (u - e ± {p) + i/(2r))~ 
rium Gq = h{uj) (G R - G$), where h{w) ■ 

Within the self-consistent Born approximation we find 
the linear (in A) correction to the Green's function, G\ . 
Having Gi we can calculate any single-particle quan- 
tity, i.e., density or current. As usual in the linear re- 
sponse theory the Keldysh component Gf- splits into 
two parts, Gf = G x ' + G x ' 1 . The first part, Gf' 1 , 
corresponds to the retarded-advanced (R-A) combina- 
tions in the Kubo formula, while G^' 11 stands for the 
R-R and A- A combinations [H|]. The standard Keldysh 
perturbation theory gives the spin-charge density matrix 
as p = ±n(q,Sl) + S (q,Sl)* = J % J [-i G<] = 

"5 / 1 1 (Cf^f - G ? + G i)> where V) is the 
charge density while s(q, f2) is the spin density. It splits 

as follows p ■ 



5' + p^where/S' = -irff J^Gf'' 



and p 



11 



(Gf < J1 -G? + G? 



1 f 2*£ f 

"2 J 271- J T2nf* [ 

Introducing the average spin-orbit band splitting 



PF\/a\ + a 2 



we can define three regimes: (i) 
"super-clean" r _1 < A 2 F m* /p 2 F = m*(a 2 i + a 2 ^); (ii) 



D' 



clean m*(a 2 ^ 



< t 1 < A F ; (iii) dirty t _1 > Aj 



Our results below apply both in the clean and in the dirty 
regimes, but not in the "super-clean" one, i.e., our results 
apply for r _1 > A 2 F m*/p 2 F . 
We obtain 



2c 



^gp B Brr 



(5) 



The functional / is defined as 
I[X{p)] = 

1 f d 2 p 



m 
1A1 



J G ° {P + 9/2 ' Ef + n/2) ' X[p) 



G${p 



(2^ 2 
q/2, E F 



n/2) , 



(6) 



while I is a 4 x 4 matrix defined by its action on the 
4-vector p as Ip — I[p] (we just use the fact that p is 
independent of p to represent the functional I[p] as a 
product of a 4 x 4 matrix / and a vector p). 

The second contribution to the density, p 11 is given by 



\gvp B B<r . Thus the total density follows from 



[1-(1 



e{vA}. 
2V~ 
iftr)!] 



2t 



gvp B Bcr 



We allow for arbitrary external frequency f2, including 



where fl > r . However, we limit ourselves to the experimen- 



tally relevant regime v F \q\ <C t _1 . Recently the spin and 
charge response functions to the longitudinal fields were 
calculated for arbitrary values of q in Ref . [l~7| • 

We expand the matrix / in powers of q, I = I^ - 1 + 
+ . . .. In zeroth order in q the matrix / is diagonal 
and its elements are given by 
1 



I, 



(o) _ 



on 



7(0) 



r(0) 
yy 



(Q + D 



r(0) = 

b 2 (l+T a 



(Q 



2aVR 
D)-b 2 {\ 



Tot) 



(8) 



2aVR 

where R = (a 2 + b 2 ) 2 - b 4 T 2 = (a 2 + (1 + T a )b 2 )(a 2 + 
(4 - T a )b 2 ), D = (a 2 + b 2 - VR) /T a , and Q = 

(a 2 + b 2 + s/ltj . We have introduced a = 4 — iClr and 
b = 2A f t. Analyzing the path of the complex function 
i?(fi) we conclude that for y/R we have to choose the 
branch-cut along the positive semi-axis R > 0. Alter- 
natively we can use the usual definition of ^/TTT (with 
the brunch-cut along the negative semi-axis) but replace 
V^R^ -iV^R. 

For the part linear in q, I^, we obtain the following 
matrix elements 

= PF|g|rsina 



zx 



,•(1) 



/(I) 



-1^ = 



mr 

PF\q\T cosa 

TO* 



(1) 

-zy 



(9) 



where 



iab 



1+T a 
R 



a 2 + b 2 (l+T a 
iab\ 



1-T Q 

R. 



a* 



b 2 (l-T a 



(40) 
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We have neglected terms mixing the charge density with 
the spin density, since they contain a small parameter 
m/(ppT) as compared to the spin-spin terms. 

Expanding the RHS of Eq. ([7]) in powers of q up to 
the terms linear in q and neglecting again the charge 
density term for the same reason as above we obtain for 
the orbital term 



iQi/I 



,{vA}, 



2c 



- r (0) rr 



r (o) 



PS,* 



+ 



ve \E\ sin a 

< 

tpf 

ve\E\\q\ m 



,(0) 



ve\E\cosa (o) 



{a)(T z , 



(11) 



where E = —A. The dimcnsionlcss cocfhcicnts in the 

c 

zcroth order in q contributions are given by 



c (0) - - - 



b(l + T a )J^(b 2 -D 



4a 



c (0) - 



6(1 - T a ) 



1+T„ 



(6 2 + D) 



4a 



(12) 



while the coefficient in the linear in q contribution is 



> z (a) = t ^T^Tl x 
(6 4 (1 - 2*) - a 4 ) 6 2 - (2T a a 2 b 4 + RD) cos 2a 



AaRVR 



(13) 



Finally, the Zeeman term in the RHS of Eq. |JJJ reads 
[1 - (1 - ifir)J] 



2r 

2t 



E [ 



1 - a/i°) 



B a a a 



(14) 



Note that the magnetic terms (oc B = x A) are already 
of first order in q. 

We are now in a position to calculate the spin density. 
First we obtain the zeroth-order in q contribution. It is 
given by 



,(o) _ 



JO) 



C 



(0) 



1 _ r (0) ^E,x 



yefi, 6(l + T ct )y / l-r ct (6 2 - D) 
2PF ' 2aVR-(Q + D) + b 2 {l+T a ) 



C 



(0) 



r(o) 



veE x 



b(i - T a )VTTT^(b 2 + D) 



JO) _ 



2pf 2aVi? - (Q - D) + & 2 (1 - T a ) ' 
0. (15) 

This is a generalization of the well known result [H, H[ 
meaning that there is an in-plane spin polarization "per- 
pendicular" to the applied electric field. We observe that, 



for T a 0, it is only perpendicular to E if the electric 
field is along one of the main axes x or y. 

Next we calculate the first order orbital contribution. 
We only calculate the z component of the spin density, 
as it was zero in the zeroth order. We obtain 



^(1), orbital 



Z r>W 

r (0) E,z 



T 


( 

1 J-zx 


1 - if) 


\~ 




[ *zy 


1 - lif 


V T 


<l + 


r(l) s (0) 
±zx &x 



Z 

r (0) ^E,x 



c 



(0) 



1 ± yy 
r(i)J0) 

J-zy By 



r(0) ^E,y 



1 

ve\E\\q\T 1 



r (0) 



m ' 



=E(«) 



i - /i? 

? -(l)J0) 



• 2 

sm a 



■(1) (0) 2 

izy'c- E ' y cos z a 



v/j, B B 2 



AO) 

±XX 

2^7 



(0) 



fm e \ l\it 
\m*J i _ /(°) 



(«) + 



•(1) (0) . 2 

izx C E x sm a 



AO) 

1 XX 



+ 



1-1, 



•(1) (0) 2 
C E, V C0S a 



AO) 



(16) 



where /ib = e/(2m e c) and m e is the bare electron mass. 
The last form is most convenient for dimensional analysis 
and plotting, as the factor in square brackets is dimen- 
sionless. 

Finally, for the Zeeman term we obtain 



(1), Zeeman 



1 - I 
1 



(o) 2r 



1 - al. 



(0) 



B z 



Ao) 

ar zz ' v g fi B 



1 



AO) 

J-ZZ 



B, 



(17) 



Discussion. We observe that both parts of the out-of- 
plane spin polarization s z , i.e., the orbital part given by 
Eq. (flBj) and the Zeeman part given by Eq. (JTTJ) can be 
regarded as linear response to the out-of-plane magnetic 
field B z . Thus our analysis amounts to a calculation of 
the susceptibility x(^,q) so that s z = x^ z - 

We observe that in the clean limit, i.e., for Ap > t -1 , 
the orbital susceptibility greatly dominates over the Zee- 
man one. As one can see in Figs. [1] and [21 for experimen- 
tally relevant parameters, the susceptibility x exceeds the 
Pauli susceptibility by a factor of order hundreds. 

In the vicinity of T a = we reproduce the results of 
Ref. [l| and the susceptibility is peaked around £1 = 2Ap. 
For T a substantially different from zero a double peak 
structure develops with the positions £1 = 2yl ± T u Af- 
Here the pole singularity present at T a = splits into 
two square-root singularities corresponding to zeros of 



4 



function R(Q) at fl = 2y / l±T a A F - i/r. One or the 
other peak are emphasized depending on the angle a as 
seen in Fig. O 



lie A 




FIG. 1: Real and imaginary parts of the total (orbital plus 
Zeeman) spin susceptibility \ = s z /B z for 2Aft — 5. Solid 
lines: T a — 0.1, dashed lines: T a = 0.5, dotted lines: 
T a — 0.9. Parameters assumed as in GaAs: m e /m* ~ 15, 
g — —0.44. The results are plotted for a — n/A which also 
corresponds to the averaged over a susceptibility. 





.A 
/ 1 

(\iA 
















FIG. 2: Real and imaginary parts of the total (orbital plus 
Zeeman) spin susceptibility \ = S z/B z for 2Aft — 5 and 
T a — 0.5. Solid lines: a = 0, dashed lines: a — ir/2, dotted 
lines: a — n/4. Parameters assumed as in GaAs: m e /m* ~ 
15, g = -0.44. 



We obtained our results for a microwave with a given 
direction of the wave- vector q, i.e., for a given angle 
a. The most obvious way to observe the orbital con- 
tribution to the spin susceptibility would be by apply- 
ing a homogeneous oscillating magnetic held B z , e.g., by 
putting the sample into a magnetic coil. Such a field cor- 
responds to an equal superposition of plane waves with 
all possible wave vectors q laying in the xy plane. To 
obtain the orbital spin response in this case one should 
just average over a, i.e., substitute (cos 2a) — and 
(cos 2 a) = (sin 2 a) = 1/2. This is what we plot in Fig.Q] 
Note, that for T a — 0, i.e., for pure Rashba or Dres- 
selhaus coupling, the response is cv-independent and the 
averaging brings nothing new. On the other hand, when 
the two couplings are of comparable strength, the suscep- 
tibility strongly depends on a (see Fig. [3]) and averaging 
over a can change the result considerably. At T a = 1 the 
orbital susceptibility vanishes. 

We thank I. Martin for numerous encouraging discus- 
sions. 
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